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EQUATION OF STATE OF DENSE MATTER AND 
MAXIMUM MASS OF NEUTRON STARS 



P. Haensel 



\£i ' Abstract. Theoretical models of the equation of state (EOS) of neutron- 

star matter (starting with the crust and ending at the densest region of 
the stellar core) are reviewed. Apart from a broad set of baryonic EOSs, 
strange quark matter, and even more exotic (abnormal and Q-matter) 
m \ EOSs are considered. Results of calculations of M max for non-rotating 

neutron stars and exotic compact stars are reviewed, with particular 
emphasis on the dependence on the dense-matter EOS. Rapid rotation 
. increases M max , and this effect is studied for both neutron stars and ex- 

m . 

of masses of neutron stars in binaries, and the consequences of such a 
confrontation and their possible impact on the theory of dense matter 
are discussed. 
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S-j 1 Introduction 

Since the pioneering paper of Oppenheimer and Volkoff (1939) the problem of 
, the actual value of the maximum mass of neutron stars, M max , is puzzling both 
J^j ' the theorists and observers. Oppenheimer and Volkoff, who used free neutron-gas 
r£ model, obtained M max = 0.72 M©. The precisely known mass of the Hulse- Taylor 
C$ ' pulsar is a direct proof of the dominating role of the nuclear forces in the equation 
of state (EOS) of neutron stars, which stiffen the EOS so that M max > 1.44 M . 
Unfortunately, the value of M max is determined by the EOS of dense matter at 
p > 2p , where po — 2.7 x 10 14 g cm -3 is normal nuclear density, and at such 
high densities the EOS is largely unknown. This results in a large uncertainty in 
theoretically derived value of M max . 

Within General Relativity, a compact object of 10 M Q > M > 1 M Q and 
R < 50 km cannot be but a neutron star or a black hole. All compact objects 
with M > -M max are necessarily black holes, while only those with M < M max 
can be neutron stars. From the observational point of view, the problem consists 
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in determining the mass function of neutron stars: its upper edge is just M max . 
However, the present measurements of neutron stars masses are most probably 
strongly biased due to a specific evolutionary formation scenario (binary radio 
pulsars) or are not very useful due to large errors (X-ray binaries; some possible 
exceptions are discussed in Sect. 1140 . Hopefully, the knowledge of neutron-star 
masses will improve in the future (see discussion in Sect. I14|) . Then, confrontation 
of measured masses with theoretical models will give us precious information on 
the EOS of matter at p ~ 10 15 g cm -3 . 

The present paper reviews the status of the theoretical determination of M max . 
Sections |2H3 ar e devoted to a brief description of existing models of the EOS of 
dense matter. Calculations of neutron-star structure in General Relativity are 
presented in Sect. El In Sects. I6~2l6. 3\ we discuss the mass-central-density diagram 
for non-rotating (conventional) neutron stars made of baryons, with particular 
emphasis on the value of Af max and its dependence on the EOS of baryonic matter, 
as well as the stability of stellar configurations. Short Sect. [7| is devoted to a 
discussion of the origin of the maximum of neutron star mass. An upper bound 
on M max is studied in Sect. |SJ The problem of stability of stellar configurations 
is reviewed in Sect. 16.31 Phase transitions in dense matter and their impact on 
neutron-star masses, and in particular on M max , are described in Sec. 03 Effect of 
rotation on neutron-star structure and on the maximum mass is discussed in Sects. 
11011 II Section El is devoted to intriguing (hypothetical) family of strange quark 
stars, and Sect. EH - to even more exotic Q-stars. Theoretical values of M max are 
confronted with existing measurements of neutron-star masses in Sect. [HI Some 
final remarks are presented in Sect. ITol 

2 Equation of state of the neutron-star crust 

We assume that the crust is built of matter in full thermodynamic equilibrium (the 
so called cold catalyzed matter corresponding to minimum energy per nucleon at a 
given nucleon density) ; the case of accreted crust will be briefly mentioned at the 
end of the present section. The crust can be divided into an outer crust, in which 
a lattice of nuclei is immersed in an electron gas, and an inner crust consisting of 
a lattice of nuclei immersed in an neutron gas and electron gas. The outer crust 
ends at the neutron drip point ~ 4 x 10 11 g cm -3 , while the inner crust extends 
down to the crust-core interface at ~ 10 14 g cm' 3 . 

Up to the neutron drip point, the EOS can be calculated using the experi- 
mental data on the neutron-rich nuclei and semiempirical nuclear mass formulae 
(Haensel and Pichon 1994). Therefore, the EOS of the outer crust is rather well 
established. As soon as one leaves the region of experimentally known nuclei, the 
EOS of cold catalyzed matter becomes uncertain. This uncertainty rises at densi- 
ties higher than the neutron drip density. The properties of nuclei are affected by 
the ambient neutron gas which contributes more and more to the total pressure. 
Therefore, the problem of correct modelling of the EOS of pure neutron gas at 
subnuclear densities becomes important. The real EOS of cold catalyzed matter 
stems from the real nucleon Hamiltonian, which is expected to describe nucleon in- 
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Fig. 1. The SLy EOS (Douchin and Haensel 2001) of neutron star crust. Dotted vertical 
lines correspond to the neutron drip and crust-core transition. 



teractions at p < 2p$. In practice, in order to make the solution of the many-body 
problem feasible, the task is reduced to finding an effective nucleon Hamiltonian, 
which would enable one to calculate reliably the EOS of cold catalyzed matter for 
10 11 g cm < p < pq, including the crust-core transition. 

We will illustrate the properties of the inner crust EOS using the results ob- 
tained by Douchin and Haensel (2001) for the SLy model of effective nuclear 
hamiltonian; it will be hereafter referred to as the SLy EOS. Notice that the 
sound velocity is v s = y/dP/dp. The overall SLy EOS of the crust, calculated 
including adjacent segments of the liquid core and the outer crust, is shown in 
Fig. ^ I n the outer crust segment, the SLy EOS cannot be visually distinguished 
from the EOS of Haensel and Pichon (1994). One notices a significant softening 
(decrease of v s ) just after the neutron drip point. At densities greater than the 
neutron drip one, the SLy EOS stiffens gradually (v s increases) with growing den- 
sity, due to the increasing contribution of dripped neutrons to the pressure. There 
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Table 1. Equations of state of the liquid core of neutron star 
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Balberg et al. 1999 



is a discontinuous increase (jump) of sound velocity at the crust-core interface. 

In the SLy EOS, the crust-core transition takes place as a very weak first- 
order phase transition, with the relative density jump of the order of a percent 
(Douchin and Haensel 2001). Let us remark that for this model the spherical 
nuclei persist to the crust bottom. For the FPS EOS, the crust-core transition 
takes place through a sequence of phase transitions with the changes of nuclear 
shapes (spheres- rods-plates-tubes-bubbles, Lorenz et al. 1993). All in all, while 
the presence of the exotic nuclear shapes is expected to have dramatic effect on 
the transport phenomena and elastic properties of matter, their effect on the EOS 
is small. 

Some neutron stars, in particular those in close binary systems, have crust 
composed of accreted matter. The EOS of accreted crust is somewhat stiffer than 
for the ground-state crust (Haensel and Zdunik 1990). However, the effect of this 
difference on the value of M max is negligibly small. 
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Fig. 2. Selected models of the EOS of neutron star cores denoted as in Table^ except for 
the BGN1 and BGN2 which here have labels BalbNl and BalbN2. Large dots correspond 
to maximum density in stable neutron stars. Asterisks correspond to the density above 
which EOS is superluminal (v a > c). Left panel: EOSs of the npep matter. Right panel: 
effect of hyperons is shown by comparing the EOS without hyperons (i.e., for the npep 
matter) and EOSs in which hyperons A, E, S are included (+H1 and +H2 correspond 
to the BGN1H1 and BGN1H2 models of Table 0. 



3 Equation of state of the neutron-star core 

At p < pq, core matter is a liquid composed mostly of neutrons with a few percent 
admixture of the equal number of protons and electrons. If the Fermi energy of 
electrons exceeds the muon rest energy (105.7 MeV), muons replace a fraction of 
electrons to minimize the energy of the system. Such a system in beta equilibrium 
is usually called the npe/it-matter. This is the simplest model of matter in neu- 
tron star cores: except for the presence of muons, which are insignificant for the 
EOS, the matter constituents - neutrons, protons, and electrons - are the same 
as in familiar terrestrial matter. Still, even for this simplest composition, the un- 
certainties in the EOS are quite large, especially at densities significantly higher 
than pq. This results from the approximations and deficiencies of the many-body 
theory of dense nucleon matter, and from the lack of knowledge of strong inter- 
actions in superdense matter. We illustrate these uncertainties in the left panel 
of Fig. |21 where four examples of the EOS of the npe/n-matter are shown. The 
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Fig. 3. Effect of phase transitions on the EOS of neutron-star core. Solid line: stable 
lower-density (1) or higher-density (h) phase. Dotted line: metastable 1-phase. The 
density at which metastable 1-phase becomes unstable with respect to the transition into 
a dense phase is denoted by p C rit. Left panel: first-order phase transition at P = Po 
from a pure 1-phase to a pure h-phase, with the density jump at P = Po from pi to ph. 
Right panel: transition from 1-phase to a mixed h-phase at P — Po and p = po. Notice 
that here po means the threshold density for the equilibrium-phase transition. The right 
panel could describe also the effect of a second-order phase transition, provided the dotted 
segment is removed. 



brief characteristics of these models and the references to the original papers are 
given in Table □ The BPAL12 and BGN2 EOSs should be considered as the 
soft and stiff extremes of the theoretical models. Only segments below large dots, 
which correspond to the maximum density in stable neutron stars calculated for 
this EOS (see Sect. 16.21 and in particular Table EJ), are relevant for neutron stars. 
Notice that the stiffest BGN2 EOS is superluminal (sound velocity v a > c) at the 
highest densities relevant for neutron stars, which reflects the inadequacy of the 
non-relativistic approach at such a high density. 

Let us denote the chemical potentials of the constituents of the npe/x-matter 
by fi n , fip, fi e , and fi^. As soon as the sum fi n + p e exceeds the in- medium energy 
of T,~ hyperon, this strange baryon will appear in as a stable constituent of dense 
matter. £~ is expected to be the first hyperon to appear in dense matter at 
some p\. The second hyperon to appear, above a higher threshold density p?, is 
A . Hyperons appear due to strangeness-changing weak interactions. At a given 
baryon density, the inclusion of hyperons lowers significantly the matter pressure 
compared to the case of the npep mater, because the highest-energy neutrons are 
then replaced by the low-energy hyperons. This effect of the softening of the EOS 
of dense matter due to the presence of hyperons is clearly seen in the right panel 
of Fig. El 

Pion or kaon condensation in dense matter, as well as deconfinement of quarks, 
predicted by some theories, imply a softening of the EOS above the phase-transition 
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threshold. The effect of such a phase transition on the EOS is visualized in Fig. El 
4 EOS of strange matter 

Strange quark matter (SQM) is a hypothetical self-bound state of deconfincd 
quarks which form a plasma of up (u), down (d), and strange (s) quarks with 
n u ~ rid — n a of total baryon number density = (n u + + n s )/S and a small 
admixture of electrons n e ~ (10~ 5 — 10~ 4 )rib. The existence of a self-bound state 
at zero pressure results from the pressure exerted by the ordinary vacuum on a 
volume of the QCD vacuum in which the quarks can move. In the Bag Model this 
pressure is equal to the bag constant B ~ 60 MeV/fm 3 . In the simplest Bag Model 
with massless free quarks the density of SQM at zero pressure is p s = AB/c 2 . The 
hypothesis that SQM is the true ground state of matter at zero pressure means 
that energy per unit baryon number -Esqm < -E( 56 Fe) = 930.4 MeV. Recent re- 
views of physics of strange quark matter are given by Glendenning (1999), Madsen 
(1999), and Weber (1999). 

As the pressure of SQM vanishes at p s (larger than normal nuclear density po), 
the EOS has a linear form P ~ [p— p s ) x constant for densities only slightly higher 
than p s . It is very fortunate that this simple linear form of the EOS turns out to be 
an excellent approximation at higher densities, up to maximum density in stable 
strange stars. This property is independent of the SQM model, and holds for the 
MIT Bag Model (Zdunik 2000) and other models of SQM (Gondek-Rosihska et al. 



Let us focus our attention on the Bag Model. A linear approximation reads 



where constants a and p s are determined by fitting the exact EOS (Zdunik 2000). 
Notice that the linear form, Eq. (|4. 1|) - is exact in the case of massless quarks (free 
or interacting). Using first law of thermodynamics one can easily see that Eq. 
l|0]| implies 



where n s is the baryon density of SQM at zero pressure. Therefore three fitting 
parameters - p s , n s , and a - completely determine an EOS of SQM. 

5 EOS of abnormal and Q-matter 

Several types of hypothetical self-bound state of dense baryon matter, which could 
constitute true ground state of the matter at P = 0, were suggested in the past. In 
1970s, it was suggested that the pion condensation could lead to the appearance 
of a self-bound, superdense state of pion-condensed nucleon matter (Migdal 1971, 
1974, Hartle et al. 1975). This phase of matter was expected to form "abnormal 
nuclei" of large A and with the density significantly higher than p (Migdal 1971, 



2000). 



P/c 2 



a(p-ps) , 



(4.1) 




(4.2) 
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1974) but no sign of experimental evidence of the existence of the "abnormal 
nuclei" has been found up to the time of this writing. 

The idea of "abnormal state" of nuclear matter, as advanced by Lee and Wick 
(1974; for a review see Lee 1975), is based on a schematic field-theoretical er-model 
of strongly interacting nucleon matter. In the abnormal state, which appears at 
sufficiently high nucleon densities, the nucleons are thought to become nearly 
massless. This is because the cr-field term couples to the nucleons as an (negative) 
addition to the nucleon rest (bare) mass, implying a nearly vanishing nucleon 
effective mass. This density-dependent effect can lead to the appearance of a 
second minimum in the n\, dependence of the energy per baryon at some n a , 
significantly higher than no, with the binding energy significantly larger than in 
the "normal state" at = n . However, the original a- model is schematic and 
does not pretend to describe quantitatively the normal nuclear matter at uq. 
Were the cr-model more complicated to give correct quantitative description of the 
nuclear matter at rib ~ no, then the abnormal state at supranuclear density would 
disappear (Pandharipandc and Smith 1975, Moszkowski and Kallman 1977). 

Supersymmetric extensions of the Standard Model of elementary particles and 
their interactions predict the existence of macroscopic self-bound superdense lumps 
of a scalar-field condensate with a well defined electric and baryonic charge - bary- 
onic Q-balls. The Q-balls were proposed as a hypothetical component of cosmolog- 
ical dark matter (Kusenko et al. 1998, Kusenko 2000). It was also proposed that 
Q-balls of a macroscopic size could have energy per unit baryon number lower than 
56 Fe. EOSs of Q-matter were constructed by Bahcall et al. (1990). A common 
feature of these models, shared with the Lee- Wick model of abnormal matter is 
that nucleons become massless inside the condensed Q-phase. Two basic parame- 
ters of the model are: the energy density Uo of the scalar field inside the Q-matter, 
and the coupling strength a v of the vector field to the nucleon. It is convenient 
to introduce a dimensionless parameter £ — a v U^ 2 Tt / \f%, and to use it instead of 
the parameter a v . 

Consider the simplest case of ( = 0. The standard value used by Bahcall 
et al. (1990) is U = 13.0 MeV fin" 3 , which corresponds to p(P = 0) = p s = 
1.0 x 10 14 g cm~ 3 : for this model the density of Q-balls is subnuclear. With 
increasing (, the EOS of the Q-matter becomes stiffer and the value of p s lower. 
In the limiting case of ( — 16 considered by Bahcall et al. (1990) (at the same value 
of U = 13.0 MeV fm~ 3 ) they get p s = 5.5 x 10 13 g cm~ 3 . The predicted density 
of self-bound Q-matter at zero pressure is two to five times lower (!) than the 
normal nuclear density, which results from a strong reduction of effective nucleon 
masses in this hypothetical state of nucleon matter. 

In all cases considered in the present section, the EOS of an abnormal phase 
of dense baryonic matter can be well approximated by 

abnormal matter, Q — matter : P ~ a(p — p s )c 2 . (5.1) 

This EOS is usually stiffer than that of strange quark matter, and one often has 
a ~ 1. 
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Finally, let us stress that there has been no observational evidence for the 
existence of hypothetical self-bound states of dense matter reviewed in the present 
section by the time of this writing. 

6 Neutron star models 

Neutron stars are relativistic objects. Their structure and evolution should be 
studied using the General Theory of Relativity. The importance of relativistic 
effects for a star of mass M and radius R is characterized by the compactness 
parameter r g /R, where r g = 2GM/c 2 — 2.95 (M/Mq) km is the gravitational 
radius. For a non-rotating black hole, r g /R = 1. Typically, one has r g /R ~ 0.2 - 
0.4 for a neutron star, while r g /R -C 1 for all other stars, e.g., r g /R ~ 10~ 4 for 
white dwarfs, and r g /R ~ 10 -6 for main-sequence stars of M ~ 1 M Q . 

We will approximate the stress tensor of neutron-star matter by that of a 
perfect fluid, i.e. non-viscous medium of total energy density £ (mass density 
p = E/c 2 ), in which all stresses are zero except for an isotropic pressure P. The 
perfect-fluid approximation is justified by the fact that the shear stresses, e.g., 
those produced by elastic strain in the solid crust or by strong magnetic field, are 
generally negligible compared to the pressure. 

The problem of finding hydrostatic equilibrium simplifies considerably in the 
case of strongly degenerate neutron star interior where thermal contributions to P 
and p can be neglected. Then, the EOS of matter involves only P and p, but not 
the temperature, a situation which is nearly always valid in neutron-star interior 
(important exceptions are: neutron-star atmospheres, newly-born neutron stars, 
envelopes of exploding X-ray bursters). 

For a static spherically symmetric neutron-star configuration the equation of 
hydrostatic equilibrium is (see, e.g., Shapiro and Teukolsky 1983) 



where m(r) is the mass contained within a sphere of (coordinate) radius r. Equa- 
tion H6.2(l describes mass balance; its apparently Newtonian form is illusive since 
the space-time is curved. Equations (|6.1|) - (|6.2fl should be supplemented by an 
EOS, P = P{p). The above equations constitute then a closed system of equa- 
tions to be solved for obtaining P(r), p(r) and m(r). 

In the stellar interior, P > and dP/dr < 0. The stellar radius is determined 
from the condition P(R) = 0. Outside the star, i.e., for r > R, we have P — and 
p = 0, and from Eq. (|6.2(1 we obtain m(r > R) = M = const. The latter quantity 
is called the total gravitational mass of the star, the total energy content of a star 
is E = Mc 2 . 

In order to construct a neutron-star model one has to solve the set of Eqs. (|6.1|) 
H6.2J1 supplemented by an EOS of neutron-star matter. Let us notice that the mass 




(6.1) 



dm 
dr 




(6.2) 
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density and baryon density are not constrained to be continuous. In principle, 
density discontinuities can appear at pressures corresponding to first-order phase 
transitions in dense matter (see Sect.[3J). The density profile is determined using 
the EOS, through p = p(P). 

Equations Htj.lfl - I|6.2|) are integrated from the stellar center with the following 
boundary conditions: -P(O) = P c , m(0) = 0. The stellar surface, r = R, is 
determined from P(R) = 0. 

EOS constitutes a crucial input for the neutron-star structure calculations. 
Models of the EOS were described in Sects. |2H3 Here we briefly describe the EOS 
used in neutron-star structure calculations discussed in the present paper. The 
outer neutron-star crust is described by the EOS of Haensel and Pichon (1994). 
The inner crust is described by the FPS model, developed by Lorenz (1991), except 
for the SLy EOS, which describes in a unified manner both the core and the crust. 
Actually, a particular choice of the crustal EOS has little importance for the global 
parameters of massive neutron-star models, which are mostly determined by the 
EOS of the liquid core at supranuclear densities. 

The EOS of the core is dominated by the contributions from strong (nuclear) 
interactions between nucleons (or, more generally, between baryons; see Sect.^J). 
The reliability of the theory of a dense, strongly interacting many-body system 
decreases rapidly with increasing p. Moreover, we have very little experimental 
information about baryonic interactions at supranuclear densities, especially about 
the interactions involving hyperons which may appear at higher densities. As a 
consequence, the EOS of neutron-star matter at supranuclear densities is strongly 
dependent on the specific microscopic theory of dense matter employed. 

Our discussion of the impact of the EOS(p > po) on the neutron-star structure 
will be based on the limited, but hopefully representative, set of models of dense 
baryonic matter. 1 The selected models of dense cold baryonic matter are listed 
in Tabled They were briefly discussed in Sect. [3] As we stressed in Sect. |3 the 
most important qualitative feature of an EOS at p > po - as far as neutron-star 
structure is concerned - is its stiffness. A simple definition of the stiffness can be 
phrased as follows: the higher the value of P at a given p, the stiffer the EOS. 
However, this definition is ambiguous and sometimes misleading. A somewhat 
better definition relies on the dimensionless adiabatic index 7 = (ny > /P)(dP/dny 1 ). 
However, a core-EOS relevant for the complete family of neutron-star models has 
to be given for po < p < (10 — 15)po- It can be rather stiff at p < (2 — 3)po, and 
then become softer at higher p, due to the hyperonization of matter with increasing 
density (see Sect. OJ- Therefore, it is useful to base the effective stiffness of the 
EOS on the value of the maximum allowable mass of neutron stars, M max , which 
is a functional of the EOS. The topmost EOS in Tabled is the softest one, and 
the effective stiffness of the EOSs increases from the top of the table to its bottom 



1 Stellar models based on the EOSs including hypothetical exotic phases, such as pion con- 
densate, kaon condensate, and deconfined quark matter, will be considered separately in Sect. 
A specific case of hypothetical strange (quark) stars or even stranger stars built of abnormal 
matter or Q-matter will be briefly reviewed in Sect. Upl and 1131 
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Fig. 4. Mass density versus radial coordinate for neutron-star models of gravitational 
mass 1.51 M© (solid line) and 1.92 M© (long dash-dotted line). Calculations are per- 
formed for the BBB2 EOS of the core, for which 1.92 M© is the maximum allowable 
mass. The positions of the crust-core interface and the neutron-drip point are indicated 
by thin vertical lines. 



(see Table EJ- Let us remind, that the stillest /softest EOSs are extreme models, 
characterized by the compression modulus of symmetric nuclear matter, which is 
significantly higher/lower than the standard "empirical value". These EOSs are 
included to represent the theoretical extremes of the EOS of dense matter. 

Two of the EOSs listed in Tabled namely SLy and FPS, are the unified ones, 
and describe, within a single physical model, both the crust and the core of a 
star (see Sect. In all other cases, an overall EOS of the stellar interior was 
constructed by matching smoothly the EOS of the crust with that of the core. 

6.1 Basic internal structure of neutron stars 

The distribution of matter within a neutron star depends on its mass and results 
from an interplay between the pressure and gravity. In Fig. 01 we show the mass 
density within a neutron star as a function of radial coordinate, calculated for the 
BBB2 EOS of the core, for two masses. The higher mass is equal to the maximum 
mass allowable for this EOS. The calculations show that the density within the 
core is rather uniform, for M = (1.2 — 1.5) M Q ; the most pronounced drop occurs 
close to the crust bottom edge. There is a steepening of the density drop near the 
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Fig. 5. Mass density versus radial coordinate for neutron-star models of gravitational 
mass 1.4 M Q , for the core described by the BPAL12 EOS (long dash-dotted line), BBB2 
EOS (solid line), and BalbN2 ( denoted as BGN2 in Table [TJ EOS (dotted line). For the 
BPAL12 and BalbN2 models, positions of the crust-core interface and the neutron-drip 
point are indicated by thin vertical lines. 



neutron drip point, which results from the softening of the EOS due to the neutron 
drip. The increase of the stellar mass increases the gravitational pull within the 
crust, which leads to the steepening of density profiles and thinning of the crust, 
most pronounced at M = M max . For a medium-stiff BBB2 EOS of the core, 
the crust contains 1.4% and 0.8% of the total gravitational mass of 1.24 M Q and 
1.51 M Q stars, and the crust thickness is 1.01 km and 0.72 km, respectively. At 
the maximum mass, the crust contains only 0.2% of the stellar mass, and the crust 
thickness reduces to about 0.29 km. The crust mass and thickness are smaller for 
softer EOSs of the core, and larger for stiffer ones. 

While the crustal EOS is rather well established (see Sect. EJ| , the crust struc- 
ture results from the interplay of its EOS and the gravitational pull exerted by the 
core: the latter depends on the core compactness. Therefore, the uncertainty of 
the EOS in the core implies uncertainty in the crust structure. This is visualized in 
Fig. [3 where we show the mass density profiles inside a 1.4 M© star calculated for 
the BPAL12, BBB2, and BGN2 EOSs of the core. The crust thickness ranges from 
0.7 km for the softest core EOS up to 1.3 km for the stiffest one. The dependence 
of the fraction of stellar mass contained in the crust on the core EOS is even more 
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Table 2. Maximum mass configurations for EOSs listed in Tabled Displayed parameters 
are: radius -R, compactness parameter r g /R, central baryon density n c , and central mass 
density p c . 



EOS 


■'"max 


R 


rjR 


n c 


Pc 




[M e ] 


[km] 




[fm- 3 ] 


[10 15 g cur 3 ] 


BPAL12 


1.46 


9.00 


0.480 


1.76 


3.94 


BGN1H1 


1.64 


9.38 


0.519 


1.60 


3.72 


BBB1 


1.79 


9.66 


0.547 


1.37 


3.09 


FPS 


1.80 


9.27 


0.572 


1.46 


3.40 


BGN2H1 


1.82 


9.53 


0.564 


1.45 


3.48 


BBB2 


1.92 


9.49 


0.596 


1.35 


3.20 


SLy 


2.05 


9.99 


0.605 


1.21 


2.86 


BGN1 


2.18 


10.9 


0.592 


1.05 


2.26 


APR 


2.21 


10.0 


0.651 


1.15 


2.73 


BGN2 


2.48 


11.7 


0.626 


0.86 


2.02 



dramatic: this fraction ranges from 0.7% for the softest BPAL12 EOS to 2.2% for 
the stiffest BGN2 one. The BBB2 model is typical for medium-stiff EOSs: the 
crust thickness is about 0.8 km and the mass fraction is 1%, respectively. 

6.2 Mass-central-density diagram and M max 

The models of cold, static neutron stars form a one-parameter family; they can 
be labelled by their central density p c , or central pressure P c . In Fig. we show 
dependence of the gravitational mass, M, on p c for some EOSs from Table ^ a t 
p c > 2.5 x 10 14 g/cm 3 . 

On the lower-density side, the M(p c ) curves exhibit a minimum M min ~ 
0.1 Mq, which depends rather weakly on the assumed EOS (Haensel et al. 2001), 
and which is not shown in Fig. 

On the higher-density side, the M(p c ) curves exhibit a maximum, which strongly 
depends on the assumed EOS, and is reached at densities a few times 10 15 g cm -3 . 
Configurations with M > M max cannot exist in hydrostatic equilibrium and must 
necessarily collapse into black holes. For the selected EOSs M max ranges from 
1.46 M Q to 2.48 M Q (Table©. 

Replacing neutron-star matter by a free (noninteracting) Fermi gas of neu- 
trons lowers the value of M max to 0.72 M (Oppenheimer and Volkoff 1939). If 
one further allows for the beta equilibrium between otherwise noninteracting nu- 
cleons, electrons, and muons, the EOS becomes even softer, with M max = 0.70 M . 
Clearly, the values of Af max obtained neglecting nuclear interactions are in con- 
tradiction with precisely measured masses of the binary radio pulsars fSect. IT4*jl . 
Turning the argument around, we may say that the measured mass 1.44 Mq of 
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Fig. 6. Gravitational mass M versus central density p c of static neutron-star models for 
different EOSs. A maximum on a M — p c curve is indicated by a filled circle. Configura- 
tions to the right of the maximum are unstable with respect to small radial perturbations, 
and are denoted by dotted lines. The shaded band corresponds to the range of precisely 
measured masses of binary radio pulsars fSect. HH . 

the Hulse- Taylor pulsar (Sect. IT4"|l implies that nucleon-nucleon interaction is suf- 
ficiently repulsive at supranuclear densities to lift M max by more than hundred 
percent from the non-interacting nucleon gas value. 

In view of uncertainties in the hyperon-nucleon (H-N) and hyperon-hyperon 
(H-H) interaction in dense matter, let us first consider a subset of EOSs with Ne/z 
composition of matter (N denotes np) . Let us remove from this subset the extreme 
models BPAL12 and BGN2. Then the subset will be restricted to the models of 
the Ne/K matter which reproduce empirical saturation properties of nuclear matter 
(see Sect.[3J). Under such conditions, we narrow the range of theoretical maximum 
allowable mass of neutron stars, 

M max (Ne/i) ~ (1.8 - 2.2) M . (6.3) 

The appearance of hyperons softens an EOS, and therefore decreases the value of 
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-Mmax- For the selected EOSs, the appearance of hyperons at p ~ 2p lowers the 
value of maximum allowable mass to an even narrower range, 

M roax (NHe//) ~ (1.5 - 1.8) M Q . (6.4) 

Let us stress that this conclusion is obtained under the assumption that the hy- 
perons appear at about 2po. However, the threshold for the hyperon appearance 
depends sensitively on the N-H interaction in dense matter, which is poorly known 
for a neutron rich dense matter at p > p . Therefore, the values of M max in Eq. 
I|6.4[l have to be taken with a grain of salt. Actually, it is reasonable to say that 
Eq. I|6.4|) corresponds to a typical effect of hyperons on M max (and on neutron- 
star structure), assuming that without hyperons M max is given by Eq. (|6.3|l . Alas, 
the lack of knowledge of the H-H interaction implies large uncertainty in the ef- 
fect of hyperons on M max . For example, an EOS calculated by Vidana et al. 
(2001) shows particularly strong softening due to the hyperonic interactions, with 
M max (Nep) = 1.89 Mq and M max (NHe/i) = 1.34 M Q . The last value contradicts 
the measured mass of the Hulse- Taylor pulsar. When Vidaha et al. (2001) re- 
move (artificially) the H-H interaction they get M max (NHe^) = 1.47 M Q , which 
is marginally acceptable. 

6.3 Stability of neutron stars 

The mass-central density diagram shows that cold, static equilibrium configura- 
tions of neutron stars can exist only for M m ; n < M < M max . The value of M m ; n 
seems to be rather well established, M m ; n ~ 0.1 M Q (see Haensel et al. 2002 and 
references therein). Equilibrium configurations of neutron stars with M < Mmj n 
do not exist. 

Non- rotating configurations with M > M max collapse into black holes. Nu- 
merical simulations of such a collapse were performed by several authors, e.g., by 
Gourgoulhon and Haensel (1993) and Baumgarte et al. (1996a, b). 

While all points on the M(p c ) curve in Fig. ^correspond to the configurations 
of hydrostatic equilibrium, only the points which correspond to stable equilibrium 
are of astrophysical interest. At a given total baryon number, hydrostatic equilib- 
rium corresponds to the extremum of M which is a functional of matter density 
distribution p; this is equivalent to the vanishing of the first-order variation of 
M, due to perturbation of p denoted as bp. However, this does not guarantee 
that the second-order variation, quadratic in dp, is positive, which is necessary 
for the hydrostatic stability. It can be shown (Harrison et al. 1965, Zeldovich 
and Novikov 1971) that stable equilibria correspond to the segments of the M(p c ) 
curve for which dM/dp c > ; this is the static stability criterion. Therefore, the 
configurations to the right of the maxima in Fig. H3 are unstable with respect to 
small radial perturbations. 2 In view of this, p c ,max is the maximum density which 



2 Actual perturbations of neutron stars are time-dependent, and can be decomposed into 
normal pulsation modes. Assuming that pulsations do not move neutron-star matter from its 
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can be reached in the interior of stable static neutron stars. The value of the 
limiting density Pcna* in static, cold neutron stars is even more uncertain than 
the value of M max : it ranges from 2 x 10 15 g cm~ 3 (about 7p ) for the stiffest EOS 
to 4 x 10 15 g cm" 3 (about I5p ) for the softest one. It should be stressed that 
even 7po is very far beyond the limits within which our models of dense neutron 
star matter can be considered as reliable and unambiguous. 



7 On the origin of M max 

When explaining the very existence of M max , two kinds of general physical argu- 
ments can be used. They are associated with expected behavior of cold matter at 
ultrahigh densities and with General Relativity. 

The arguments based on the high-density behavior of the EOS of dense matter 
were proposed, for Newtonian stars, by Landau (1932). His brief paper presented 
a general physical argument why above some M max a large self-gravitating sphere 
of cold matter cannot sustain itself against gravitational collapse. 3 Harrison et al. 
(1965) rephrased the arguments in a way applicable to neutron stars. However, it 
should be stressed that this way of reasoning is based on the assumption that in the 
high-density region in which the stability is lost the EOS of dense matter is that 
of the ultra-relativistic free Fermi gas. The counterexample to such a behaviour 
was presented by Zeldovich (1962). 

Irrespectively of the form of the EOS of dense matter, the upper bound on M 
is a consequence of the General Relativity. Consider the right-hand-side of the 
OV equation of hydrostatic equilibrium, Eq. (|6.1I) . It describes gravitational pull 
acting on a unit proper-volume of matter, 

, „ Gmp ( 4ttP\ ( P\( 2Gm\~ 1 , N 
gravitational pull = ^- 1 H T 1 + — =■ 1 5- . 7.1 

The gravitational pull is given by a Newtonian-like term —Gmp/r 2 , multiplied by 
three relativistic factors. With increasing M, all three factors in brackets amplify 
the pull compared to the Newtonian case, and the increase of M with an increase 
of central pressure P c becomes harder. Mathematically, the derivative dA//dP c 
decreases with growing M. 

Let us illustrate this property with an unphysical case of incompressible fluid of 
mass density p nlc = constant. The total gravitational mass is then (47r/3)/9i„ c i? 3 , 
and the pressure profile within the star is determined analytically (solution was 
obtained by Karl Schwarzschild in 1916, and is described in Box 23.2 of Misner et 



ground state (i.e., that the EOS for both static and pulsating configurations is the same) one 
can show that the configurations to the right of M max in Fig.|H]are indeed unstable with respect 
to the fundamental mode of small radial pulsations (Harrison et al. 1966). The finite timescale 
of reactions between constituents of neutron star matter which may be longer than the pulsation 
period complicates the dynamical analysis of stability of neutron stars with respect to small 
radial pulsations (Meltzner & Thorne 1967, Chanmugan 1976, Gourgoulhon et al. 1995). 
3 Landau's paper was submitted for publication before discovery of neutron. 
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al. 1973). The central pressure P c can be determined as a function of r g /R. As 
the mass increases, P c as well as r g /R increase monotonically, too. For P c — > oo, 
the radius tends to a finite value i?M max = | r g , an d therefore the mass tends to 

= = ^ 1/2 = 5.09 M f 5xl ° 14gCm ' 3 ) ' 

9 G 35/ 27r l/2 G 3/2 p l/2 ° V p inc J 

(7.2) 

This limit is the General Relativity effect: there is no limit to the mass of the 
incompressible-fluid stars in Newtonian gravitation. If Af max exists for an incom- 
pressible fluid, then a fortiori it should exist for any EOS of dense matter with 
finite compressibility. However, the numerical value of M max is mainly determined 
by the EOS for p > 2po which is largely unknown. 

8 Upper bound on M max 

In view of the uncertainty in the value of M max , it is of interest to have a possibly 
firm upper bound on M max . Let us assume that the EOS is known up to a certain 
density p m . This reliably known segment of the EOS will be denoted by P < (p). 
It spans the range < P < P u = P < (p u ). Having fixed P<(p), we can treat 
M max as a functional of the EOS at P > P u ; this unknown EOS will be denoted 
as P > (p). The inverse function p > (P) 1 determined for P > P u , does not need to 
be continuous. Let us require only that the P>{p) is subluminal and fulfills the 
condition of stability: 

0<, s 2 = ^<c 2 , (8.1) 
dp 

where v s is the speed of sound in the dense medium. Our task is to find the 
maximum of the functional M max [P>(p)] on a set {P>(p)} of EOSs which satisfy 
conditions l|8.1[l . As it turns out, M max is maximized by the so called causal-limit 
(CL) EOS, 

P^(p) = P u + (p- Pu )c 2 , (8.2) 

which therefore yields the upper bound we are looking for, M^ ax . It is quite 
natural: the CL EOS is the stiffest one at p > p u . Numerical calculations show 
that for p u < 2p the effect of the outer layers with p < p u on the value of M^ x 
is negligibly small. Actually, had we used the pure causal-limit EOS of the form 
P = (p — p s )c 2 (where p s is the surface density of a star built of matter with such 
an EOS), we would slightly increase (by less than one percent) the value of M max 
compared to that obtained for the CL EOS with p u = p s . 4 This can be easily 
explained, because P<(p) is not maximally stiff. All in all, for p u < 2po one gets 
(Rhoades and Ruffini 1974, Hartle 1978, Kalogera and Baym 1996, Glendenning 



4 For a pure CL EOS one gets an exact result M max = 3.0 (5 X 10 14 g cm- 3 /p s )V 2 M Q . 
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1999) 

1 

Vs <c =^ M max <M^ = 3.o( 5Xl ° ^ gCm " V M . (8.3) 

The maximum allowable mass for any EOS(p > p u ) with subluminal sound velocity 
is lower than the above upper bound. As the inequality l|8.3J) is widely accepted, it 
seems to be safe to say that actual M max of static neutron stars built of baryonic 
matter is below 3 M . 5 A rapid rotation increases the value of M^ x , as discussed 
in Sect. El 

9 Phase transitions in dense matter, third family of compact stars, and 
maximum mass 

The existence of exotic phases of dense matter in neutron star cores is a subject of 
debate. Unfortunately, since the physics of dense matter at p > 2po is uncertain, 
the problem cannot be solved within the present many-body theory. Therefore 
one is forced to admit that the question: is a specific exotic phase of dense matter 
present in a neutron star core? can only be answered by unambiguous identifica- 
tion of signatures of this phase in neutron star observations. 

In what follows, we will distinguish the case of a moderate and a strong soft- 
ening of an EOS by a phase transition. We will refer to a phase transition as 
moderate, if it does not lead to the appearance of an unstable segment in the 
M — p c curve. In the opposite case of a strong softening of the EOS, the phase 
transition will imply the existence of an unstable segment of the M — p c curve, 
which will separate a family of lower-density configurations from a distinct fam- 
ily of superdense compact objects. Let us start with the case of the moderate 
softening of the EOS, without density jump (right panel of Fig.|3J). Then p(P) is 
continuous, A = Ph_/ p\ = 1, but the adiabatic index 7 = (ny,/ P)dP/dnb suffers a 
jump (a drop) at the appearance of a new phase at P = Pq. The M — p c curve is 
displayed in Fig. E^,. Such a situation corresponds to hyperonization, and also to 
transitions to a mixed phase (baryons and pion or kaon condensate, baryons and 
dcconfincd quark matter). In the latter case, the mixed phase exists up to some 
limiting density at which a transition from a mixed phase to a pure higher-density 
phase takes place; it is is accompanied by a second discontinuous change (actually, 
an increase) in 7. 

Consider now a phase transition which implies a strong softening of the EOS 
without a density jump. This is the extreme case of a very strong decrease of 7 
after the phase transition discussed above, followed by a hardening of the EOS at 
still higher densities. In this extreme situation the M — p c curve has two stable 



5 Notice that we assume that neutron star has an outer layer of density p < p u , composed 
of normal cold dense matter. For some very exotic hypothetical models of compact stars with 
superdense surface built, e.g., of a self-bound Q-matter, the limit of 3 Mq does not apply (Sect. 
1131 . 
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Fig. 7. Mass - central-density relations for an EOS, for which a phase transition at a 
pressure Po implies softening with no density jump (right panel of Fig. 3). Configura- 
tion Co is the last one composed exclusively of a lower-density phase, its central pressure 
P c = Po- Thin lines represent configurations calculated using the EOS without a phase 
transition. Thick lines show configurations calculated for the EOS with a phase transi- 
tion. Dotted segments correspond to configurations which are unstable with respect to 
small radial perturbations. Left panel (a): moderate softening of the EOS. Right panel 
(b): strong softening of the EOS, implying appearance of an unstable branch between 
Ci^ and C^in (thick dotted line), and a separate branch of superdense stars between C^m 



branches of static equilibrium configurations, Fig. 03, corresponding to two distinct 
families of neutron stars. The first (lower-density or normal) family is continuously 
connected with low-mass neutron stars and terminates at C* r The second family 
contains "compact neutron stars", and will be referred to as the "higher-density 
family" or "superdense family" . The configurations belonging to the superdense 
family have central densities p^ ain < p* < p* nax and masses M^ in < M* < M^ ax . 
They are more compact and tightly bound than those containing the same number 
of baryons and belonging to the normal branch. 

The appearance of a distinct family of superdense stars deserves a general 
remark. Superdense stars form actually the third family of compact stars built 
of degenerate matter, the remaining two being the well known white dwarfs and 
"lower-density" neutron stars. 

Two examples of the strong softening of the EOS without density jumps are 
given by Glendenning and Kettner (2000) and Schacffncr-Biclich et al. (2002). In 
the model of Glendenning and Kettner (2000) the softening of the EOS results from 
the appearance of a quark-baryon mixed phase and takes place for a substantial 
quark fraction. At higher densities, their EOS stiffens when the mixed phase is 
replaced by a pure quark phase. In the model of Schaeffner-Bielich et al. (2000) 
the softening is due to a copious appearance of hyperons. It is obtained for a 
specific model of hyperon-hyperon interaction. It is followed by a substantial 
stiffening of the EOS at higher densities. Similarly to the first example, it leads to 
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Fig. 8. Mass - central-pressure relations for an EOS containing a phase transition with a 
density jump (left panel of Fig. 3) . Configuration Co is the last one composed exclusively 
of a lower-density phase, with P c = Pq. Notations are the same as in Fig. Q Panel 
(a): moderate softening of the EOS. Panel (b): strong softening of the EOS, implying 
appearance of an unstable branch between Co and C^ in (thick dotted line) , and a separate 
branch of superdense stars between C^ in and Cj^ ax . For further explanation see the text. 



the appearance of the second stable "higher-density branch" of the M — p c curve; 
the models belonging to this branch are mostly composed of hyperons. In both 
examples, the maximum mass of the "higher-density branch" is lower than that 
of the "lower-density" family, M* lax < M max . 

At the same baryon number, the energy M*c 2 of the high-density configura- 
tion is lower. Let A^ in and A^ ax denote the limiting baryon numbers of sta- 
ble configurations on the higher-density branch. Configurations with extremal 
mass and extremal baryon number coincide (see, e.g., Harrison et al. 1965). For 
A^ in < A < A^ ax , the true stable configurations are those belonging to the higher 
density branch, while the configurations consisting of the same number of baryons 
and lying on the lower-density branch are metastable. 

As a third type of phase transition in dense matter, we consider a transition 
with a moderate softening associated with a density jump in the EOS. This jump 
occurs at a pressure Pq , at which a pure lower-density 1-phase of density p\ coexists 
with a pure higher-density h-phase of density ph, so that the relative density jump 
is A = ph/p\ > 1 (left panel of Fig.OJ. In view of the discontinuity of p c at Pq, the 
equilibrium configurations have to be parametrized by the central pressure P c . 

One has to stress a particular role played by the quantity 

^, s f(i + A). m 

A moderate softening associated with a first order phase transition corresponds to 

A < A cr it- 

Finally, we will consider the strong softening of the EOS associated with a 
first order phase transition with A = ph/pi > A cr it . In this case, the appearance 
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of a small core of h-phase destabilizes the neutron star. Specifically, for static 
configurations C* with a small h-core and A > A crit one has dM/dp c < 0. Such 
configurations are therefore unstable with respect to small radial perturbations and 
collapse into stable configurations with a large h-core. The instability condition A > 
|(1 + Po/pic 2 ) was first derived by Seidov (1971), using the static energy method. 
Ten years later this condition was rediscovered by Kaempfer (1981) who studied 
the necessary condition for the onset of the neutron-star collapse initiated by the 
phase transition at its center. It is worth to mention that the Newtonian version of 
this criterion (A > §) was first derived by Lighthill (1950, see also Ramsey 1950) 
in the context of stability of two-phase planets. The relativistic effects stabilize 
neutron stars with small h-cores by AA cr i t = A cr it — § = |Po/pic 2 , which can be 
as high as ~ 0.2. Static stable equilibrium configurations of neutron stars split 
into two families, visualized in Fig. |Ht>- The superdense branch C^C^ax I0rms 
a third family of compact stars, apart from the white dwarfs and lower-density 
neutron stars. It should be stressed, that in contrast to the instability implied 
by a second-order phase transition due to the hyperonization or appearance of a 
mixed quark-baryon phase (Sect.yJ, a typical situation for A > A cr it corresponds 
to M max < M^ ax and A max < A^ ax (see, e.g., Brown and Weise 1976, Haensel 
and Proszyhski 1982, and Migdal et al. 1990). 

10 Rotating neutron stars 

We limit ourselves to a stationary rigid rotation. In this case the angular frequency 
Q of rotation of a matter element, measured by an observer at infinity, is constant. 
A rotating configuration has axial symmetry with respect to the rotation axis. 
A stationary rotation of stellar bodies in General Relativity has been studied by 
many authors; the present state of this field is reviewed in a Living Review by 
Stergioulas (2001). 

Rotating configurations form a two-parameter family, and can be labeled, e.g., 
by the values of p c and SI. These configurations are flattened and their equatorial 
radius, R cq , is larger than the polar radius, i? po i. Configurations C(p Cl fl) cover a 
region in the M — R cq plane. A curve M(i2 eq , f2) for a fixed f2 is limited on the 
high- and low-density sides. On the high-density side, the curves are limited by 
the condition of stability with respect to small axi-symmetric perturbations. In 
order to check whether a configuration C(p c i,J7i) with angular momentum J\ is 
stable, one has to consider a family of configurations C(p c , J\) with fixed angular 
momentum J = J\ in the neighborhood of C{p c \, J7i). Configuration C(p c i, fii) is 
stable if 



where the derivative is calculated along the {C(p c , J±)} family. A line determined 
by (dM/dp c )j = separates the configurations stable with respect to small ax- 




(10.1) 
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isymmetric perturbations from the unstable ones. 6 

The low-density boundary is determined by the condition of stability with 
respect to the mass-shedding from the equator. A necessary condition for the 
existence of a stationary rotating configuration is the following: the equatorial 
velocity of an element of stellar matter has to be smaller than the velocity of 
a test particle moving on a circular orbit of radius i? eq in the equatorial plane, 
the so-called Keplerian velocity Uk- The Keplerian velocity corresponds to the 
Keplerian angular frequency fix = U-^/ R cv called also the mass- shedding angular 
frequency. 

In order to illustrate the effect of rotation on the structure of observed pulsars, 
consider the radio pulsar PSR B1937+21 with shortest observed period effect of 
rotation i^g = 1.56 ms (Backer et al. 1982); the corresponding angular frequency 
is fi^ax — 641 Hz. The effect of rotation on the maximum-mass configuration 
is very small. For M ~ M max the period of 1.56 ms corresponds to the slow 
rotation regime (fi -C fii<), an d therefore Mj^ ms — is quadratic in the 

small parameter fi = 0/ y/GM/R 3 . Consider for example the specific case of 

the SLy EOS (Douchin and Haensel 2001). Then ff = (fi^ ax ) 2 - 0-06. As will 
be shown in Sect. 1111 the highest rotation frequency allowed by the condition of 
stationarity implies the increase of M max by some 20%. Therefore, the fractional 

obs 

increase of M max connected with rotation at P = 1.56 ms is 0.2(fi max ) 2 ~ 2%, 
which agrees very well with the exact numerical result (Haensel & Douchin 2001). 
The effect will be smaller for a softer EOS (e.g., FPS EOS) and larger for a stiffer 
EOS (e.g., APR EOS). 

11 Maximum mass of rotating neutron stars 

Rotation increases M max because the centrifugal forces oppose the gravity. For 
fi <C fix, this increase is quadratic in fi/fij< and therefore very small. More 
generally, M max (fi) — M^ ax is an even function of fi, because it does not depend 
on the orientation of the spin axis. 7 

Within the entire set of stable stationary configurations {C(p c , fi)} one can find 
two important extremal configurations: with the maximum mass M™ ax and with 
the minimum period P m i n . These two configurations do not coincide but are very 
close to each other. Depending on the EOS, the maximally rotating configuration 
with the rotation period P m in can have the central density higher or lower than the 
maximum-mass configuration, but the difference is at most a few percent (Cook 
et al. 1994). As a rule, the mass of the maximally rotating configuration is lower 
than M™ ax by less than one percent (Cook et al. 1994). 



6 In the static limit, J\ = 0, we recover the condition of stability with respect to radial 
perturbations (which are a special case of axisymmetric perturbations), dM/dp c > (see Sect. 

7 We remind that our stars are built of an ideal fluid and the effects of the magnetic fields are 
neglected. 
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It is useful to note that for realistic baryonic, subluminal (v s < c) EOSs, 
M™^ is (to a very good approximation, typically within 3%) proportional to the 
maximum mass of non-rotating configurations (Lasota et al. 1996) 

AC* X ~ 1.18J<£ . (11.1) 

However, the above formula is not valid for the EOSs of the form P/c 2 = a(p— p s ), 
and in particular, for strange quark stars; their case will be considered in Sect. 1121 
In Sect. we derived an absolute upper bound on the static neutron-star mass 
based on the knowledge of the EOS at p < p u ~ 2po under the constraint of v s < c. 
Let M^^ stat be the upper bound for non-rotating neutron stars. Rotation will 
increase the upper bound, M^ K (fl) > M^f tat . The upper bound for rotating 
stars, M^^ rot , is obtained for the same causality-limit EOS as for the non-rotating 
models; it is reached at ^ very close to O max . Its precise value was obtained by 
Koranda et al. (1997) 

vs<c: Mff. < M%£* = 3.89 M 2 , (11.2) 

For a given p u , it is 30% larger than the static upper bound, M^£ tat . 

The rapidly rotating configurations, considered in Sects. fTOVll II are stable with 
respect to the axisymmetric perturbations and mass shedding. However, they can 
be susceptible to various secular instabilities, reviewed by Stergioulas (2001). 



12 Strange quark stars 

By strange stars we will mean the compact objects built entirely, or predominantly, 
of self-bound strange quark matter (SQM, Sect.0}. The possibility of the existence 
of self-bound strange quark stars built entirely of SQM, was contemplated by 
Witten (1984), who considered a simplified model of SQM, composed of massless 
u, d, and s quarks, with the EOS of the form P = \{pc? — B) (see also Brecher 
and Caporaso 1976). Witten showed that for the bag constant B sa 60 MeV fm -3 
close to the value needed to reproduce experimental masses of baryons within the 
MIT Bag Model, the parameters of the maximum-mass configuration for strange 
stars are similar to those for realistic neutron stars built of baryonic matter. The 
first detailed models of strange stars, based on a more realistic EOS of SQM, 
taking into account strange quark mass and the lowest-order QCD interactions, 
were constructed by Haensel et al. (1986) and Alcock et al. (1986). These authors 
discussed the basic physical properties of strange stars, and their astrophysical 
manifestations. 

Further development of physics and astrophysics of strange stars focused on the 
refinement of the EOS of SQM (particularly, beyond the MIT Bag Model), and on 
specific properties of strange stars, such as neutrino emissivity, rotational proper- 
ties, superfluidity, pulsations, electromagnetic radiation, and cooling. Physics and 
astrophysics of strange stars is reviewed by Glendenning (1999), Weber (1999), 
and Madsen (1999). 
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The surface density of bare strange stars is equal to that of SQM at zero 
pressure, p s . It is therefore some fourteen orders of magnitude larger than the 
surface density of normal neutron stars. 

In what follows, we will illustrate the generic properties of strange stars assum- 
ing SQM1 EOS of SQM of Zdunik et al. (2001). The M-R curve for bare strange 
stars is shown in Fig. ED For strange stars with M > 1 M Q the radius changes very 
little with M, R ~ 9—11 km. It is quite similar to that obtained for neutron stars 
with a moderately stiff EOS. However, at lower masses, the radius of bare strange 
stars behaves in a completely different way. Namely, it decreases monotonically 
with decreasing M, and R oc Ms for M < 0.3 M . Such a behavior of low-mass 
bare strange stars can easily be explained within the bag model. Gravitational 
pull decreases rapidly with decreasing M, and can be neglected at M < 0.3 M 
compared to the pressure of the normal vacuum on the volume "filled" by the QCD 
vacuum: this pressure confines SQM to a sphere of radius R. Due to very high 
incompressibility of strange matter, the density within a low-mass strange star is 
nearly constant and close to p s (see Fig. ITU)) . On the other hand, the low-mass 
strange stars can be described in the Newtonian theory, which gives M ~ ^-p s R? 
and Roc Mi. 

The internal structure of bare strange stars is very different from that of neu- 
tron stars. First of all, their surface density is huge, p s ~ 10 15 g cm~ 3 . The 
density profile in such a star is very flat (Tig. I10|). Even at the maximum mass, 
i.e., under the maximum gravitational compression, the central density is only five 
times higher than the surface one; this is to be compared with fourteen orders 
of magnitude center-to-surface density difference for neutron stars! The density 
difference decreases rapidly with decreasing M. For a 1.4 M Q bare strange star the 
central density is only forty percent higher than the surface one. For low-mass bare 
strange stars of M < 0.3 M , the center-to-surface density difference is negligibly 
small. 

12.1 Maximum mass of strange quark stars 

As mentioned before, the EOS of SQM is accurately determined by two param- 
eters, a and p s (Sect. 0J. Consider bare strange quark stars. Using the linear 
representation of the EOS and passing to dimensionless variables P, p, r, 



one can rewrite the equations of hydrostatic equilibrium of strange stars in a di- 
mensionless form, At a fixed a, the equilibrium configuration can be obtained from 
a universal dimensionless solution of transformed Eqs. i|tj.l[) - (|6.2|) by returning to 
ordinary variables r, P, p. This property of Eqs. H6.1f) - I|6.2fl implies useful scaling 
properties relating a solution which corresponds to an EOS P = a(p — p' s )c 2 to 
the solution for the EOS P — a(p — p s )c 2 ■ In particular, the maximum mass for 
different values p s and p' s are related by the scaling relation 




,2 ' 



~ P 
P = — , 
Ps 




(12.1) 
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Fig. 9. Mass-radius relation for bare strange stars and for strange stars with normal crust 
calculated for the SQM1 EOS. Extremal (maximum and minimum mass) configurations 
are indicated by filled circles. Strange stars with crust have maximum crust mass and 
thickness: density at the crust bottom is equal to the neutron-drip one, 4.3 x 10 11 g cm -3 . 
The dotted segment represents configurations unstable with respect to small radial per- 
turbations. The hatched area corresponds to the region of the M — R plane prohibited 
by General Relativity and by condition v s < c. Long dash-dot line gives the radius of the 
marginally stable circular orbit around a strange star (it is the radius of the innermost 
stable orbit of a particle orbiting in the equatorial plane of a non-rotating star). The 
vicinity of the minimum-mass configuration for strange stars with the crust is shown in 
the right panel. 



M^ ax = M max • j , (12.2) 
For a strange matter EOS based on the MIT Bag Model, the scaling relation, 
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Fig. 10. Mass density p versus radial coordinate r for three bare strange stars of different 
masses, calculated for the SQMf EOS of strange quark matter. The highest of three 
masses is the maximum allowable mass for this EOS. 



Eq. I|12.2[l . turns out to be precise within better than one percent (Haensel et al. 
1986). Similar precision is obtained for other models of SQM (Gondek-Rosinska 
et al. 2001). 

Scaling properties become particularly simple (and actually exact) for the MIT 
Bag Model EOS with non-interacting massless quarks, i.e., for a = | and p s = 
AB/c 2 . This EOS will be referred to as SQMO. The formula for M max for the 
SQMO EOS was derived by Witten (1984), 

1 96 

M max = -== Mq , (12.3) 

V -°60 

where B ea = B/(60 MeV fm" 3 ). 

The scaling formula enables one to explain why the parameters of massive 
strange stars (M ^1 — 2 Mq) are so similar to what we expect for normal neutron 
stars. Consider first the simplest SQMO model. For the hypothesis of strange 
matter to be correct, energy per unit baryon number in SQM should be less than 
930.4 MeV which is the minimum energy per nucleon in normal baryonic matter 
at zero pressure and temperature; this value is reached for the 56 Fe crystal. This 
condition implies B 60 < 1.525. On the other hand, a consistent model of SQM 
should not lead to a spontaneous fusion of neutrons into droplets of u, d quarks, 
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which would eventually transform into droplets of SQM - "strangelets" . This 
in turn means that the energy per unit baryon number in the u, d matter with 
2n d = n u should be above 939.6 MeV, which implies B 60 > 0.982. All in all, 
both constraints lead to 1.6 M < M max < 2.0 M Q , which reminds us a typical 
maximum mass range calculated using medium-stiff EOSs of baryon matter. An 
allowance for the finite mass of the strange quark and for the QCD interactions 
within the MIT Bag Model can increase maximum masses of strange stars due to 
a simultaneous decrease of the allowed values of B (this results from the condition 
E S qm < 930.4 MeV, see Zdunik et al. 2000). 

A strange star can be covered by a solid crust of matter with density lower than 
the neutron drip point (outer neutron-star crust). The nuclei of the crust do not 
fuse with the SQM core because they are separated from the SQM by a repulsive 
Coulomb barrier. However, since the mass of the crust is lower than 10 -4 M Q , its 
effect on M max of strange stars is negligible. 

12.2 Maximum mass of rotating strange quark stars 

Rotation increases the maximum mass of strange quark stars stronger than Af max 
of neutron stars. This difference results from the different matter distribution 
within strange quark stars: the density profile is relatively flat and the surface 
density is huge. This feature increases the effect of centrifugal forces on the stellar 
structure. For simplest SQM0 EOS one gets the exact result (Gourgoulhon et al. 
1999) 

SQM0 EOS : M™* x = 1.44M£j£ . (12.4) 

In the case of more realistic EOSs of SQM relation i|12.4fl holds only approxima- 
tively. Still, on may say that the increase of M max for strange quark stars due to 
rapid rotation (about 40%) is twice that for neutron stars (20%). 

13 Maximum mass of Q-stars 

The models of hypothetical Q-stars were constructed by Bahcall et al. (1990). As 
mentioned in Sect. [SJ two basic parameters of the model are: the energy density Uq 
of the scalar field inside the Q-matter and the coupling strength a v of the vector 
field to nucleon. The results for Q-star parameters show a simple scaling with Uq 
provided the dimensionless parameter £ = a„[/g 2 tt/\/3 (Sect.jSJ is kept constant. 
Then, the scaling of M max with Uq is the same as that with respect to p s for strange 
quark stars, Eq. I|12.3I) . but with different numerical coefficients. For < C < 16 
and Uq = 13.0 MeV/fm 3 , one gets an astonishingly high maximum allowable mass, 
4.0 < M max < 8.3 Mq. This stems from a low value of Q-matter density p s at zero 
pressure: remember that nucleons are nearly massless in the Q-matter. Consider 
the simplest case of £ = 0. Then, the EOS for the Q-matter coincides with that 
for strange matter of massless, non-interacting quarks. However, the standard 
value used by Bahcall et al. (1990) while constructing the families of the Q-stars 
is Uq — 13.0 MeV fm -3 , which corresponds to p s = 1.0 x 10 14 g cm -3 , only one 
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third of the normal nuclear density! The maximum allowable mass can then be 
calculated using Eq. i|12.3fl . and is 4.0 M Q , in agreement with Figs. 4 and 8 of 
Bahcall et al. (1990). With increasing £, the EOS of the Q-matter becomes stiffer. 
In the limiting case ( = 16 considered by Bahcall et al. (1990) (at the same value 
of U = 13.0 MeV far 3 ) they get p s = 5.5 x 10 13 g cm" 3 and Af max = 8.2 M . 
This is not surprising: for a = 1 (v s = c) we can use a pure-causal limit formula 
of Sect. |H1 getting an estimate of the maximum mass « 9.0 M Q . This estimate 
reproduces within 8% the exact value of M max , obtained by Bahcall et al. (1990). 
Rapid rotation could further increase the maximum mass of Q-stars by some 30%, 
to well above 10 M©. 

Summarizing, while the arguments for the existence of Q-matter stem from 
sophisticated supersymmetric extensions of field-theoretic models of dense nucleon 
matter, the practical reasons for shockingly high Ajf max for Q-stars are very simple. 
Namely, the predicted density of self-bound Q-matter at zero pressure is three to 
five times lower (!) than the normal nuclear density, which results from a strong 
reduction of effective nucleon masses in this hypothetical state of matter. 

Additional complication is that the typical average density of Q-stars is sig- 
nificantly lower than the nuclear density. Therefore, a conversion of a neutron 
star into a Q-star should be accompanied by a significant inflation of stellar size, 
which can be obtained only at the expense of a gigantic work done against the 
gravitational pull. This should be a very peculiar type of a transformation, in 
which a compact star becomes less dense, but more bound, because nearly all 
rest mass of nucleons annihilated in the phase transition. Finally, while the exis- 
tence of supersymmetric (Q) ground-state of matter may be not in conflict with 
the terrestrial nuclear physics, reaching this state during stellar evolution may be 
virtually impossible due to huge energy barrier separating the Q-matter from the 
normal state of dense matter. 

14 Confronting theory with observations 

Some neutron stars belong to binary stellar systems. In these cases, one can try 
to measure their masses by analyzing orbital motion of the binary. Up to now 
more than one hundred such binary systems containing neutron stars have been 
discovered which are quite different in nature and observational appearance. In 
what follows we briefly review evaluations of neutron star masses based on the 
analysis of binaries containing neutron stars. 

The accuracy of measuring masses of X-ray pulsars or X-ray bursters is rather 
poor because of many obstacles in obtaining high-precision X-ray data, difficulties 
in establishing the parameters of the orbital motions in the binary and accounting 
for many interfering factors (e.g., accretion flows, tidal forces). Direct mass esti- 
mates have been obtained for several X-ray pulsars and are displayed in Fig. HP 
They will be hereafter denoted as Mj. 

In what follows, we focus our attention at the binaries with highest Mx- Recent 
analysis of the mass of neutron star in an X-ray binary, Cyg X-2 (which is an X- 
ray burster), is - at the la (68%) confidence level M x {la) = 1.78 ± 0.23 M Q . 
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Simultaneously determined mass of the companion star in Cyg X-2 is M c = 0.60 ± 
0.13 M Q (Orosz and Kuulkers 1999). However, in order to use the value of M x 
to constrain dense matter EOS, it seems reasonable to require at least 2a (i.e., 
95%) confidence level, which corresponds to M x (2cr) = 1.78 ± 0.46 M Q . This 
would generate 1.32 M Q as a lower limit on M max , which is not really useful for 
constraining modern EOSs of dense matter. With additional constraint on the 
companion mass, resulting from theoretical models (M c > 0.75 M Q ), Casares et 
al. (1998) get M x > 1.88 M Q . This would be very restrictive as far as the EOS is 
concerned, but, in view of additional assumptions and strong model dependence, 
it cannot be used as a clean measurement of neutron star mass. 

Another X-ray binary, studied for a long time, is Vela X-l (X-ray pulsar). The 
central value of Mx obtained by different authors is high, ~ 2 M© , but the errors 
are unfortunately quite large. In view of this, the lower bound on the mass of 
Vela X-l is never significantly higher than 1.4 M Q . Very recently, Barziv et al. 

(2001) obtained M x {2a) > 1.54 M Q and M x {3a) > 1.43 M Q . The 3cr lower 
bound has (formally) similar significance as the binary - radio-pulsar masses, but 
is unfortunately close to the mass of the Hulse- Taylor pulsar, and therefore does 
not yield a new constraint. 

Other examples of determination of neutron-star masses of the X-ray pulsars 
can be found in the paper of van Kerkwijk et al. (1995). One has to conclude, that 
lower limits on M max , determined at the 2cr confidence level from the condition 
M m&x > Mx , are at present too low to yield useful constraint on the modern dense 
matter EOS. This situation may change in the future. Very recently, Clark et al. 

(2002) analysed the high-mass X-ray binary 4U1700-37 (in which the compact star 
is neither an X-ray pulsar nor an X-ray burster) and found Mx — 2.44 ±0.27 M , 
and Mx (2cr) > 2.0 M Q . If the compact object is a neutron star then its mass 
would rule out soft and moderately stiff EOSs (in particular EOSs which include 
hyperons and phase transitions). 

Observations of kHz quasi-periodic oscillations (QPOs) in low-mass X-ray bina- 
ries could (potentially) yield interesting constraints on neutron-star masses in these 
binary systems, if this phenomenon is connected with the orbital motion of matter 
around neutron star. However, as by the time of this writing (July 2002), the very 
mechanism of the phenomenon of the kHz QPOs is not established. Therefore, 
resulting constraints on neutron-star masses are ambiguous. 

Fortunately, there are close binaries containing a radio pulsar and a compact 
companion (another neutron star or a white dwarf). The presence of a radio 
pulsar has a great advantage over purely X-ray sources since radio observations 
enable one to measure the pulsar periods P with extremely high precision (more 
than 12 correct digits!). Pulsars are known to be excellent timers: their proper 
rotational periods (in comoving reference frames), in many cases, are wonderfully 
stable. If, however, a pulsar participates in orbital motion, its period, as measured 
by a distant observer, varies due to the Dopplcr effect. Actually, the radio pulses 
observed on Earth are also influenced by the spacetime curvature in the vicinity 
of the companion star. Measurable relativistic effects include periastron advance, 
gravitational redshift and transverse Dopplcr shift in the orbit, and orbital decay 
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PSR J0437-4715 
PSR J1012+5307 
PSR J1045-4509 
PSR J1713+07 
PSR B1802-07 
PSR J1B04-2718 
PSR B1855+09 
PSR J2019+2425 



PSR B1518+49 

PSR B1518+49 companion 

PSR B1534+12 

PSR B1534+12 companion 

PSR B1913+16 

PSR B1913+16 companion 

PSR B2127+11C 

PSR B2127+11C companion 

PSR B2303+46 

PSR B2303+46 companion 



PSR J0045-7319 



1 2 

Neutron star mass (M G ) 



Fig. 11. Neutron star masses from observations of close binaries containing radio pulsars. 
Error bars indicate 68% (la) confidence limits. Upper limits, indicated by arrows, arc 
one-sided 95% (2a) confidence limits. Estimates of masses of eight neutron stars in four 
double neutron star systems are shown at the top of the diagram. In two cases, the 
average neutron star masses in a system is known with much better accuracy than the 
individual masses; these average masses are indicated with open circles. Estimates of nine 
neutron stars in neutron star-white dwarf binaries are shown in the central and bottom 
part of the diagram. The estimate of of the neutron-star mass in a neutron star-main- 
sequence star binary is shown at the bottom of the diagram. Vertical lines are drawn 
at M = 1.35 ± 0.04 Mq, and delimit the 68% bounds resulting from the the maximum 
likehood Gaussian distribution of the measured neutron-star masses, with mean mass 
1.35 M and standard deviation 0.04 M Q . From Thorsett and Chakrabarty (1999). 



due to the emission of the gravitational radiation. Under favorable geometrical 
conditions, Shapiro delay and geodetic precession can also be measured. A review 
of the present status of the timing analysis of the binary pulsars is given in a Living 
Review by Lorimer (2001). 

The masses of six neutron stars in the three neutron-star binaries contain- 
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ing radio pulsars have been determined with high accuracy. In the case of PSR 
B1534+12, PSR B1913+16 and their neutron-star companions, the masses are 
known within 0.2% and 0.02%, respectively. Unfortunately, such double neutron- 
star systems are expected to be formed in a very specific evolutionary scenario, 
and therefore the most precisely measured neutron-star masses most probably do 
not give much information about the complete neutron-star mass function. 

An EOS of dense matter has to explain the measured neutron-star masses. A 
simple rule is to exclude those EOSs which predict lower M max than the highest 
precisely measured neutron star mass (in July 2002: 1.44 M Q ; in our opinion result 
of Clark et al. (2002) has to be carefully checked as far as its model-dependence is 
concerned before being used to rule out soft and moderately stiff EOSs) At present, 
this condition excludes only the softest EOSs appearing in literature (some EOSs 
of dense matter with hyperons: see, e.g., Pandharipande 1971, Balberg and Gal 
1997, Vidana et al. 2000). All of recently developed EOS considered in the present 
paper are consistent with the Hulse-Talor pulsar mass (see Fig. EJ). 

Needless to say, it would be highly desirable to accurately measure the masses 
of more massive neutron stars, the more massive the better. The highest measured 
neutron star mass -M^™ implies the observational constraint 

M max (EOS) > M^ s ax) , (14.1) 

where M max (EOS) is maximum allowable mass for static neutron star models for 
an EOS. Rotation at periods longer than say 5 ms has almost no effect on M max . 
A definite discovery of (1.8 — 2.0) Mq neutron star would rule out EOSs with 
hyperons and/or phase transitions; the candidates for such massive neutron stars 
are Vela X-l and compact star in 4U1700-37, discussed previously in this section. 
A discovery of, say, 2.1 Mq neutron star would leave us with very stiff EOSs of 
dense matter containing nucleons only. 

15 Concluding remarks 

A theorist who strongly believes in the power of the theory could find a loophole 
in the criterion expressed by Eq. (|14.1|) . Why not imagine two families of "neu- 
tron stars" , described by different EOSs which correspond to two different forms 
(phases) of dense matter. Let the first of them has maximum mass Mm2x and 

(2) 

the second one A/max- So one can can contemplate a situation, in which the first 
EOS violates criterion 114. lfl but is not eliminated by it because the most massive 
observed "neutron star" belongs to a second family (Haensel and Zdunik 1989). 
For example, one might in principle consider the existence of strange quark stars 
and of neutron stars, or even of Q-stars and neutron stars. Of course, a necessary 
condition for some contact with reality is the stability of both these families. This 
means, that at the same baryon numbers compact star configurations belonging 
to different families have to be separated by a sufficiently high energy barrier. 

How far can one go with theoretical considerations not substantiated by solid 
experimental basis? Many theorists (to which the author belongs) think that re- 
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specting the Occam's razor principle is necessary in research work. 8 The fact that 
a model is not ruled out by observations is not a proof of its reality. The model 
has to be necessary for understanding the observations to be considered as repre- 
senting a reality. It would be intellectually arrogant to believe that the universe 
is filled with objects predicted by the theories based on a distant extrapolation 
of the laboratory physics. Conventional models of neutron stars should be con- 
sidered as long as they are sufficient to understand astrophysical phenomena and 
measurements. If one day a compact object of, say, 8 M Q is discovered which is 
not a black hole, then it cannot be but a Q-star. Future observations will hopefully 
bring us more information, and enrich our knowledge of the compact stars, but for 
the time being the principle devised by William of Occam is a useful complement 
to the theory of dense matter. 
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